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Abstract 

The differences between uni-directional and bi- 
directional polymerization are considered. The uni- 
directional case is discussed in the framework of the 
RNA world. Similar to earlier models of this type, 
where polymerization was assumed to proceed in a 
bi-directional fashion (presumed to be relevant to peptide 
nucleic acids), left-handed and right-handed monomers 
are produced via an autocatalysis from an achiral sub- 
strate. The details of the bifurcation from a racemic 
solution to a homochiral state of either handedness 
is shown to be remarkably independent of whether 
the polymerization in uni-directional or bi-directional. 
Slightly larger differences are seen when dissociation is 
allowed and the dissociation fragments are being recycled 
into the achiral substrate. Key Words: RNA and DNA 
polymerization, enantiomeric cross-inhibition, origin of 
homochirality. 

Introduction 

The origin of homochirality is usually believed to be 
closely connected with the origin of life (see Bada 1995 
for an overview). It may have even been a prerequisite for 
life in that the structural stability provided by chiral poly- 
mers may have been essential for the assembly of the first 
replicating molecule. If this is so, it would probably mean 
that the origin of homochirality had to be a physical one. 
Possible candidates for a physical origin of homochirality 
include the presence of polarized light from a nearby neu- 



tron star (Rubenstein et al. 1983), magnetic fields (Thie- 
mann 1984, Rikken and Raupach 2000), or mechanisms 
involving the electroweak force (e.g., Hegstrom, 1984). 
However, Bailey et al. (1998) and Bailey (2001) showed 
later that supernova remnants have not actually displayed 
circularly polarized light. Another perhaps more likely 
possibility is that homochirality developed rather as a con- 
sequence of life. This would mean that some primitive 
form of life should have been possible without chirality 
having played any role in this. 

In connection with the origin of life one used to dis- 
cuss the hypothesis of a relatively simple self-replicating 
molecule (e.g. Frank 1953). This picture ignores the pos- 
sible importance of compartmentalization that may be re- 
quired for achieving the concentrations necessary for the 
chemical reactions to take place. This led to the concept 
of a very early lipid world that would have preceded the 
often discussed RNA world. Some insight into these ideas 
can be gained by looking at recent theoretical attempts 
to build life from scratch invoking a series of steps and 
chemical processes that are thermodynamically possible 
(Rasmussen et al, 2003). Interestingly enough, their ap- 
proach involves peptide nucleic acid (PNA) because of its 
charge carrying properties and the molecule's hydropho- 
bic backbone. Its potential as contemporary genome, 
which would for example require a machinery for pro- 
tein transcriptase, was not utilized at this stage, although 
it may undoubtedly become a candidate for carrying ge- 
netic information at later evolutionary stages. 

Although this is speculation and details are unknown, 
the idea of a combined PNA/lipid world provides an at- 
tractive scenario for discussing the origin of homochi- 



1 



rality in the context of genetic evolution (Nelson et al, 
2000, Pooga et al, 2001). We picture here a situation 
where PNA has developed to having autocatalytic prop- 
erties, just like RNA in the RNA world (Woese 1967). 
PNA can be achiral if its peptide backbone is derived from 
glycine. The step toward a chiral backbone invoking, for 
example, poly-alanine, seems like a relatively minor one 
(one of two H is being substituted by CH 3 in the CH 2 
piece). However, there are two different ways of doing 
this, leading to either L or D alanine. The assembly of 
mixed L and D PNA poly-nucleotides is unlikely, just as 
it is unlikely in the corresponding case of DNA polymer- 
ization (Joyce et al, 1984). Moreover, the addition of a 
nucleotide of opposite handedness is known to 'spoil' fur- 
ther polymerization (also known as 'enantiomeric cross- 
inhibition'). This makes it increasingly unlikely to gen- 
erate L and D polymers of any appreciable length greater 
than just a few. 

The main difference between PNA and DNA polymer- 
ization is that DNA can attach new monomers only on 
the 3' end of the ribose sugar (e.g. Turner et al, 2000), 
so polymerization is uni-directional and can only proceed 
in one direction. By contrast, PNA does not have this 
restriction and can polymerize in a bi-directional fashion 
in either direction. The latter case has been addressed in 
a number of recent studies starting with Sandars (2003), 
but the former case is more readily amenable to laboratory 
verification, as is shown by recent experiments confirm- 
ing the process of enantiomeric cross-inhibition (Schmidt 
et al. (1997) and Kozlov et al. 1998). Given that the differ- 
ences between uni-directional and bi-directional polymer- 
ization have not yet been explored, we must first extend 
the formalism of Sandars (2003) to the uni-directional 
case and focus then on the comparison between the two. 

Although enantiomeric cross-inhibition seems to be an 
important ingredient of homochirality, this can only work 
if the production of new monomers is somehow biased to- 
ward the enantiomeric excess of the already existing poly- 
mers - even if this bias is extremely tiny. This is the sec- 
ond important ingredient of homochirality and is known 
as autocatalysis. Certain chemical reactions are indeed 
known to have such properties (Soai et al, 1995, Sato et 
al. 2003, Mathew et al, 2004). It is important to point 
out that these reactions are only based on dimerization, 
but they are nevertheless quite valuable in establishing the 
basic elements of homochiralization in chemical systems 



(Kitamura et al. 1998, Plasson et al. 2004), and can lead 
to quantitative predictions (Blackmond et al. 2001, Buono 
and Blackmond 2003). For a recent review see the paper 
by Mislow (2003). 

The importance of the combined action of enantiomeric 
cross-inhibition on the one hand and autocatalysis on the 
other has been well known since the very early work of 
Strong (1898) and a seminal paper of Frank (1953), who 
first proposed a simple mathematical model consisting of 
only two variables representing the relative numbers of 
left and right handed building blocks. His paper was 
tremendously insightful in that he understood not only 
the two basic ingredients needed for homochirality, but he 
was also aware that there are two rather different scenarios 
through which homochirality can be achieved, depending 
basically on how frequent the creation of a potential life 
bearing molecule is. If the creation of life was sufficiently 
frequent, life may have emerged at different locations on 
the Earth's surface (including the oceans), giving rise to 
the interesting possibility of having different life forms 
of opposite handedness simultaneously. This is the case 
studied recently by Brandenburg and Multamaki (2004), 
who estimated that left and right handed life forms could 
have coexisted for not more than the first 500 Million 
years. This is because different populations will spread 
over the Earth's surface and come eventually into contact, 
extinguishing one of the two life forms. The other possi- 
bility is that the creation of life was an infrequent event, in 
which case there was ever only one life form, which was 
then the one that led eventually to the global population 
over the Earth's surface. Regardless of which of the two 
scenarios applies, the final outcome would have been the 
same. 

In his paper, Frank (1953) only analyzed the second 
alternative in detail. This is also the scenario discussed 
in most of the approaches since then, which all discuss 
homochirality as the result of a bifurcation process [see 
also Saito and Hyuga (2004a) for a recent classification of 
different possibilities]. This forms also the basis for the 
model discussed in the present paper, where we present 
a modification of a detailed polymerization model pro- 
posed recently by Sandars (2003). In this model the enan- 
tiomeric excess grows exponentially in time. However, 
if the creation of life is a frequent event, the process to- 
ward global homochirality can only occur linearly in time 
(Brandenburg and Multamaki 2004; see also Saito and 
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Hyuga 2004b for related work). 

In the model by Sandars (2003), autocatalysis is incor- 
porated by assuming that the rate of monomer produc- 
tion of given handedness is proportional to the concen- 
tration of polymers of the same handedness. As noted 
above, this effect alone, i.e. without the additional effect 
of enantiomeric cross-inhibition, cannot lead to complete 
homochirality, because the initial enantiomeric excess is 
not (or only weakly) amplified. In order to model this 
quantitatively, Sandars (2003) assumed that polymeriza- 
tion can, at a certain rate, also occur with monomers of 
opposite handedness. This reaction produces chemically 
inactive products and it acts thus as a means of removing 
oppositely oriented building blocks (that are already in the 
minority) from the system. This model has been studied 
further by Wattis and Coveney (2005) and by Branden- 
burg et al. (2005a, hereafter referred to as BAHN) who 
showed that, for large enough fidelity of the catalyst, the 
departure from a homochiral state occurs exponentially 
fast at a growth rate that depends on the fidelity and the 
rate of enantiomeric cross-inhibition. They also discussed 
a model consisting only of primers and dimers which can 
be reduced to a set of two ordinary differential equations 
which are similar to those of Frank (1953). An impor- 
tant difference to Frank's model is the form of the cross- 
inhibition term. As discussed by Blackmond (2004), the 
feedback term in his model corresponds to the forma- 
tion of inactive dimers with one left and one right handed 
building block. This is unrealistic, because dimers with 
two left or two right handed building blocks should also 
form. This let her to the conclusion that the dimers must 
act as catalysts. 

We have emphasized that the original model of San- 
dars assumed that polymerization can occur on either end 
of the polymer. While this may be a reasonable assump- 
tion in general (and probably also for PNA), it is not re- 
alistic for RNA polymerization where polymerization can 
usually only proceed in a uni-directional fashion. Since 
uni-directional polymerization leads to a simpler model, 
and since the derivation of the bi-directional polymeriza- 
tion model has already been discussed elsewhere (see, 
e.g., BAHN), the uni-directional case is ideal for intro- 
ducing the basic ingredients of the model. Following the 
mathematical description of the uni-directional model, we 
present numerical solutions that show that the main con- 
clusions obtained from the earlier bi-directional polymer- 



ization models carry over to the uni-directional case. This 
addresses possible objections that the Sandars model is 
not applicable to RNA and DNA polymerization that is 
more easily amenable to detailed laboratory verification. 

Polymerization model 

The starting point of the model is a basic polymerization 
process 

L n + L\ — L n+ i, (1) 

where L n denotes left handed polymers of length n and 
ks the reaction rate. The corresponding model of the 
polymerization process reads 

^W = -fe[£i] ([£„]-£„_!), (2) 
at 

where [L n ] is the concentration of L n . New building 
blocks are continuously added to the model, e.g. by the 
inclusion of a substrate that provides a source Ql of new 
monomers, i.e. 

d N 

-[L 1 ] = Q L -J2k S [Li}[L n ] (3) 

n— 1 

The solution of Eqs. (|2} and (|3} is simply a wave 
traveling toward longer polymers at velocity ks[L{\ (see 
Fig.[0, as can also be seen by considering the continuous 
limit of this equation, d[L n ]/dt = — ks[Li]d[L n ]/dn. 
Note that, in contrast to a similar result for bi-directional 
polymerization (see BAHN, their Fig. 1), the functional 
form of [L n ] is continuous between n = 1 and n = 2. In 
the bi-directional case [L{\ is about twice as large as [L2]. 

The model becomes more interesting when the right 
handed polymers, R n , are also included. The interaction 
between the mirrored strands is assumed to occur through 
two separate phenomena: enantiomeric cross-inhibition 
and enzymatic autocatalysis. The autocatalysis makes the 
left handed, respective right handed, polymers catalyze 
the production of left, respective right, handed building 
blocks. The source terms Ql and Qr are proportional to 
the concentration of the achiral substrate [S] and a corre- 
sponding reaction coefficient kc- In the case of perfect 
fidelity, / = 1, the source terms are written as 

Q L = k c [S}C L , Qn = kc[S]C R , (4) 
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Figure 1: Wave-like propagation of a finite amplitude 
perturbation in the uni-directional polymerization model. 
The initial profile is a gaussian. Note the undisturbed out- 
ward propagation of the wave at n = N. The time dif- 
ference between the different curves is 20/(fc s Q) 1 / 2 . We 
have shown the first and last times as dashed and solid 
lines, respectively, and all other times as dotted lines. The 
parameters are N = 50 and kc/ks = 1. 

where Cl and Cr are some measures of the catalytic ef- 
fect of the already existing left-handed and right-handed 
polymers. The should be a monotonically increasing 
function of the overall concentration of the left handed 
polymers. The exact functional form of these expressions 
are not known. In fact, different authors have chosen dif- 
ferent prescriptions for Cl and Cr. The qualitative re- 
sults of the models do however not seem not affected by 
this choice. We find it natural to assume that 



JV 



Cl = ^2 n[L n ], 



N 

Cr = n \ R - 
n=l 



(5) 



(6) 



In the more general case of finite fidelity of the assumed 
autocatalysis, i.e. for < / <, we model there will be 
'cross-talk' between between the two handednesses, so 
we write 



Ql = k c [S] { \ (1 + f)C L + \ (1 - f)C R + Col } , 



Q R = k c [S}{l(l + f)C R + ±(l- f)C L + Cor), (8) 

Here the terms Col and Cor allow for the possibil- 
ity of non-catalytic production of left and right handed 
monomers. However, in the following we assume Col = 
Cqr, = 0. (The inclusion of Col and Cor terms leads to 
so-called imperfect bifurcations; see Fig. 6 of BAHN.) 

The enantiomeric cross-inhibition occurs when a build- 
ing block attaches to a polymer of the opposite handed- 
ness. The resulting polymer cannot continue to grow at 
the affected end and can therefore be considered spoiled. 
This phenomenon has been observed in experiments by 
(Joyce et ai, 1984) who studied template-directed poly- 
merization. When the cross-inhibition is included, the set 
of reactions in the model is (for n > 2) 



Li, 



Rx 



fes 
ki 



L n Ri, 



(9) 
(10) 



and for all four equations we have the complementary re- 
actions obtained by exchanging L and R. The new pa- 
rameter kj measures the rate at which the cross-inhibition 
occurs. The rate equations now read (for n > 2) 



d[L n ] 
dt 

d[R n ] 
dt 



k s [Li] ([£ n _i] - [Ln]) - fc/[L„][fli], (11) 



= k S [Rl] {[Rn-l] - [Rn]) - fc/[i?n][il]- (12) 



(7) 



The evolution of the spoiled polymers, L n Ri and R n Li, 
can be discarded, because, in contrast to bi-directional 
polymerization, their concentrations do not enter the uni- 
directional model. 

In comparison with bi-directional polymerization we 
note that here for n = 2 there is no extra 1/2 factor in front 
of the [Li] 2 and [Ri] 2 terms in Eqs. (O and O- This 
is because with polymerization from either end the total 
reaction rate would be twice as big. However, when two 
monomers interact, the corresponding reaction equation 
is the same for uni-directional and bi-directional poly- 
merization, because the two reacting monomers are in- 
distinguishable. Thus, whether the first binds to the sec- 
ond or the second to the first monomer does not make a 
difference. This is then equivalent to saying that for two 
monomers polymerization can occur both on the 3' and on 
the 5' end of the ribose sugar. In effect, this removes an 
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awkward 1/2 factor for the n = 2 equations in the model 
of Sandars (2003); see also Eq. (7) of BAHN. 

The reactions (|9jl and dlOt imply the presence of addi- 
tional loss terms in the evolution equations of monomers, 
so instead of Eq. (|3} we now have 



dt 

d_ 

dt 



[Li] = Q L -X L [L 1 ], 
[Ri] = Qr-\r[Ri], 



(13) 
(14) 



where we have defined decay rates 



N 



JY 



X L = ks[ [Lx] + ^[L„] +k r J2[Rn], (15) 



n=l 



N 



N 



(16) 



Comparing again with the bi-directional model, the 
present model has an extra [L{\ (or [Ri]) term, but there is 
no factor of 2 in front of the ks and kj terms and the sums 
over the concentrations of semi-spoiled polymers are also 
absent. 

From symmetry considerations it follows that there al- 
ways exist a racemic steady state ([R n ] — [L n ]) of the 
rate equations. In fact, we can show that a steady state is 
given by (for n > 2) 



-(n-1) 



[Li] (racemic), (17) 



In particular, if ki = ks, then [L n ] = [Li]/2 n ~ 1 , i.e. [L n ] 
drops by a factor of 2 from one n to the next. This is also 
true between [Li] and [Ri], while in the bi-directional 
model their ratio is 4. 

While the existence of a racemic solution is trivial, 
the interesting question is whether there exist other fixed 
points of the equations, and in this case which of these 
fixed points are stable under certain conditions. As was 
shown in BAHN the model typically goes through a 
pitchfork bifurcation from a single stable fixed point (the 
racemic solution) to a state with two homochiral stable 
fixed points where the racemic solution corresponds to an 
unstable fixed point. The order parameter controlling the 



bifurcation is the fidelity / of the autocatalysis. In Fig.|2] 
we show the enantiomeric excess, defined here as 



V 



Cr — Cl 
Cr + Cl 



(18) 



for ki/ks = 1 and kj /ks = 0.1. We also compare with 
the corresponding result from the bi-directional polymer- 
ization model. The difference between the two cases is 
however surprisingly small. 

Polymer dissociation 

The model described in the last section provides a possi- 
ble explanation of homochirality, without appealing to ex- 
ternal mechanisms for the symmetry breaking. One may 
also argue that the model is rather realistic in that it ex- 
plicitly considers the polymerization process. Less satis- 
factory are some of the details in the description of the 
polymerization process. Perhaps most importantly, the 
polymerization process is irreversible, no chain-breaking 
is included in the model. As we have already pointed 
out in an earlier paper (Brandenburg et al, 2005b, here- 
after referred to as BAN), this is unrealistic because for 
large enough fidelity the polymer length always tends to 
diverge. Also, the model cannot be self-contained since 
there is no feedback from the polymers back to the sub- 
strate. 

Before discussing in more detail the differences be- 
tween uni-directional and bi-directional polymerization in 
the presence of dissociation, let us first recall the main as- 
pects of the polymerization model with dissociation, as 
developed recently by BAN. The dissociation process is 
described by the reaction 



L, 



is 



Lr, 



and the corresponding reaction for the right handed poly- 
mers. It turns out that there are a number of subtleties that 
need consideration when constructing the detailed model 
of the chain breaking. For example, if we assume that 
the fragments can continue to polymerize, the result is a 
catastrophic over-abundance of the short chains. The rea- 
son for this is that all building blocks (L\ and R\) are used 
to produce longer polymers whereas polymers of length 
two or more cannot (according to the reactions above) 
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Figure 2: Bifurcation diagram for two different values 
of ki/ks (=1 in the upper panel and 0.1 in the lower 
panel). Note the transition from a racemic to homochi- 
ral state as a function of the autocatalytic fidelity /. Ho- 
mochirality is measured in terms of r\ = C£2 n n({L n ] — 
[Rn])/(J2 n n ([L n ] + [R n })- For weak enantiomeric cross- 
inhibition (kj/ks = 0.1 in the lower panel) the range 
of permissible values of the fidelity parameter is de- 
creased, demonstrating the importance of enantiomeric 
cross-inhibition. 



agglomerate into longer polymers. One way of remedy 
would of course be to include the agglomeration in the 
model, but the disadvantage of this is that the model then 
becomes significantly more complex due to the higher de- 
gree of nonlinearity. These issues are discussed in fur- 
ther detail in BAN, where also a number of possible of 
the model model are considered. We focus here on the 



Figure 3: Isotactic equilibrium states with polymeriza- 
tion, dissociation, and recycling of fragments into the 
substrate, for different values of M (upper panel), and 
the mean polymer length Nl (lower panel, solid line), 
compared with the bi-directional polymerization model of 
BAHN (dotted line). 



model where the polymerization fragments are recycled 
back into the achiral substrate. In the rest of this paper we 
discuss the modifications necessary to incorporate disso- 
ciation in a uni-directional polymerization model. 

In the presence of dissociation, the new system of equa- 
tions is 



(n - 1)75 [L n ], 
p^-(n-l) ls [R n ], 
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where p^ and indicate the terms due to polymer- 
ization described above. The source term in the substrate 
is given by 

Q = W L + W R + W LR + W RL (19) 

where 

N N 

W L = Y. UW n L) > W R = nW n R) ' (20) 
n— 1 n— 1 

is the total number of recycled building blocks (both left- 
handed and right-handed), and 

N N 

W LR = £(n + 1)«4 £ *> , W RL = ^(n+ l)wi R V 

n—1 n—1 

(21) 

are the corresponding contributions from fragmented (in- 
active) polymers. 

Like in the bi-directional case, the average polymer 
length scales with a quarter power of the parameter 
M = {kshs) E„=i n ([£„] + [Rn])- Thus, in order to 
achieve appreciable polymer length, the normalized total 
mass must be sufficiently large. 

Histograms of the chain distribution and the depen- 
dence of the chain length on the total normalized mass 
are given in Fig. [3] and compared with the bi-directional 
case. For small chain mass (M < 10) the chains tend 
to be very short (JVx, « 1...2), which is common to both 
bi-directional and uni-directional cases. For larger total 
mass, however, the two cases begin to depart from each 
other such that for the same total mass the chains are 
slightly shorter in the uni-directional case. 

Conclusions 

In the present paper we have modified the polymerization 
model of Sandars (2003) such that polymerization is only 
possible on one of the two ends of the polymer. Although 
PNA polymerization is probably still bi-directional, this is 
normally not the case for RNA polymerization. The sig- 
nificance of considering RNA polymerization is that it is 
readily amenable to direct experimental verification (e.g., 
Joyce 1984). One of the perhaps most curious properties 
of the model is the wave-like evolution of the polymer 



length after initializing the polymerization process. This 
prediction could possibly be tested experimentally by set- 
ting up a range of different polymerization experiments 
that are being stopped at different times. A subsequent 
analysis, as it is done for DNA sequencing, might then 
reveal a structure as seen in Fig.^ 

We emphasize that homochirality appears sponta- 
neously when two separate mechanisms are present in the 
polymerization process: autocatalysis and enantiomeric 
cross-inhibition. The accuracy of the autocatalysis is pa- 
rameterized by a fidelity factor. At low fidelity the poly- 
merization leads to a racemic solution whereas at higher 
fidelity a homochiral state is reached from an initially (al- 
most) racemic solution. The corresponding bifurcation 
diagram displays a classic pitchfork bifurcation and the 
autocatalytic fidelity acts as a control parameter. The dif- 
ferences between uni-directional and bi-directional poly- 
merization are however surprisingly small. 

In the second part of this paper we have extended 
the model to include dissociation within the framework 
of uni-directional polymerization. As in the case of bi- 
directional polymerization, the model becomes chemi- 
cally more realistic in that longer chains are now possible. 
Moreover, the model is constructed to be self-contained in 
that the need for external replenishing of the substrate is 
now replaced by the recycling of dissociation fragments. 
With respect to chirality, the qualitative behavior of the 
model is shown to persist the inclusion of dissociation. 
We therefore conclude that the existence of a transition 
between a racemic and homochiral state, as a function of 
the autocatalytic fidelity, is a robust phenomenon within 
the class of models under consideration. 
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